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Abstract: Four-wave mixing near resonance in an atomic vapor can pro- 
duce relative intensity squeezed light suitable for precision measurements 
beyond the shot-noise limit. We develop an analytic distributed gain/loss 
model to describe the competition of mixing and absorption through the 
non-linear medium. Using a novel matrix calculus, we present closed-form 
expressions for the degree of relative intensity squeezing produced by this 
system. We use these theoretical results to analyze experimentally measured 
squeezing from a ^^Rb vapor and demonstrate the analytic model's utility 
as an experimental diagnostic. 
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1. Introduction 



Relative intensity squeezing by four-wave mixing in an atomic vapor is emerging as a promis- 
ing technique for performing high-precision measurements beyond the shot-noise Hmit. First 
demonstrated by McCormick et al. [1], the technique uses atomic coherences to produce quan- 
tum correlated "twin beams", enabling the shot-noise of one beam to be measured and sub- 
tracted from the other to obtain a low-noise differential measurement; for example of a weakly 
absorbing sample. This scheme was recently shown to reduce the relative intensity noise by 
9.2±0.5dB below the shot-noise limit [2], and noise reduction has been observed in both the 
low Fourier frequency [3] and multi-mode imaging [4,5] domains. 

Furthermore, as one of the twin beams is near-resonant with the atoms, this squeezing tech- 
nique has promising applications in quantum information processing [6,7]. However, absorp- 
tion near resonance degrades the quantum correlations. Both mixing gain and absorption losses 
occur simultaneously as the beams propagate through the vapor, and are therefore competing 
processes. 

Earlier theoretical investigations of this system have applied numerical methods [3] and the 
Heisenberg-Langevin formalism [8] to predict the resulting degree of squeezing. The numerical 
model demonstrated excellent agreement with experimental results, but it can be difficult to gain 
insight into the competing processes from numerical calculations. The Heisenberg-Langevin 
model provided a microscopic description of a specific four-wave mixing configuration in a 
cold atomic gas, which accurately predicted the resulting gain profiles. However, calculation of 
the predicted squeezing required complex matrix integrals and no comparison to experimentally 
measured squeezing was presented. 

In this work, we present a very general approach for determining the squeezing produced by 
a four-wave mixing system, and develop a matrix-based analysis method to include arbitrarily 
many injected vacuum modes. Considering special cases, simple closed-form expressions are 
easily obtained. Finally, we present experimentally measured squeezing from four-wave mixing 
in a rubidium-85 vapor, and demonstrate how the model can be used as a diagnostic tool to 
determine the limiting technical factors. 

2. Relative intensity squeezing 

The "double-A" four-wave mixing scheme introduced by McCormick et al [1] uses a high- 
intensity "pump" beam to drive a cycle of four off-resonant transitions in a hot rubidium va- 
por, causing the emission of correlated "probe" and "conjugate" photons (Fig. 1 A). The probe 
transition is stimulated by a seed laser incident at an angle 9 to the pump, resulting in the 
spontaneous emission of the conjugate on the opposite side of the pump beam (Fig. IB). The 
beam powers are measured individually and subtracted to obtain the relative intensity noise as 
measured on a spectrum analyzer (S.A.). 

Labelling the Fock-space annihilation operators of the probe, conjugate and pump by a, b 
and c respectively and the interaction strength by ^ , the interaction picture Hamiltonian is 

In the "undepleted pump" approximation, the intense pump beam remains in its initial coherent 
state I \j/c) and the substitution c —?' y/c can be made: 

Mi=iTi{^ yr^.Pa' -^*{\ii* fab). 

The time-evolution of this Hamiltonian over the interaction time-scale T is 



S = exp(— /J^T//i) = exp(i^^fl^ — s*ab), where s = <^</c T. 



(1) 




This is the two-mode squeezing operator for modes a and b, where s is the "squeezing parame- 
ter" [9]. The four-wave mixing system therefore produces a two-mode squeezed state, reducing 
amplitude difference noise at the expense of increasing phase difference noise [5]. 

The phase of s results in a rotation of the (arbitrary) measurement quadratures, so s may be 
taken as real and positive. The probe and conjugate modes a and b are then transformed as 

fl 5^fl5 = cosh(i)fl + sinh(i)£^ and P ^ b^ S ^ ?,m\\{s)a + cqs\\{s)1)' . (2) 

Defining the number operator of the incident probe beam as .^Vq = ajflo and making the bright 
beam approximation (.^Vq) ^ 1, the number operators after squeezing become 

(M,> = (fl^fl) ^ G{No) and (A^^) = {Pb) ~ (G - 1 ) {No) , 

where G = cosh^ s is the increase in probe intensity, termed the "mixing gain". 

The relative intensity operator Na — Ni, is unchanged by S, so Var (Na — N/,) = Var (.^Vq) . 
Hence the beams have been amplified without increasing the relative intensity noise; they are 
relative intensity squeezed. The noise figure of the process (or "degree of squeezing") is the 
ratio of the measured noise to the corresponding shot-noise level for equal optical power. As- 
suming an initially shot-noise limited probe, the noise figure is 

Var(M,-A^.) _ {No) 

{Na) + {Nb) G{No) + {G-l){No) 2G-\ 

Therefore the measured noise power can be reduced arbitrarily below the shot-noise limit in 
the limit of ideal detection. However, optical losses are unavoidable and occur both within the 
medium (e.g. absorption) and after it (e.g. imperfect detection). These losses randomly eject 
photons from the probe and conjugate beams, decorrelating their intensities and degrading the 
observed degree of squeezing. We now construct models to quantify this effect. 

3. Optical losses after squeezing 

We initially consider only losses that occur after mixing, such as from imperfect optical trans- 
mission or detection efficiency. These losses are modelled by a beamsplitter with an empty 
port [10] whose output state is a combination of the input and vacuum modes, contributing 
"vacuum fluctuations" to the transmitted beam [11]. Denoting the vacuum modes introduced 
by losses on the probe and conjugate by the annihilation operators x and y respectively, the 
standard beam-splitter input-output relations [12] give 

^,a+^J\-riaX and b ^ ^,b + - 7]by, (4) 



where rja and rji, are the fractions of the probe and conjugate intensities transmitted. The relative 
intensity noise can then be expressed in terms of the individual beam variances and covariance 
to give 

Var(M, -A^,,) = T7,?Var(M,) + I7«(l - 17„)(^«) + rilWai{Nb) + I7fo(l - rih){Nh) 
-2riarihCoVai{Na,Nb) . 

Computing the variances using Eq. (2), the noise figure corresponding to four-wave mixing 
followed by optical losses is 



NF-ll 2(G-l)(G(r?«-T7,)^-%^) 
Gri, + iG~l)rib 



(5) 



This expression highlights the importance of balanced beam detection, as unbalanced losses 
(rja 7^ rji) result in detection of amplified noise instead of squeezing. 

4. Optical losses during squeezing: Interleaved gain/loss model 

The four-wave mixing process consists of Raman transitions between the hyperfine ground 
states (Fig. lA), which are most efficient when the intermediate virtual level is tuned close to 
resonance. However, this also increases direct absorption from the Doppler broadened transi- 
tion, increasing losses and reducing correlations. To analyze this trade-off, we develop a model 
for the effect of competing mixing and absorption on relative intensity squeezing. 

Following the approach of the numerical model presented in Ref. [3], the beam trajectories 
through the medium are divided into discrete interleaved stages of gain and loss (Fig. 2). 
Distributed models of this type were first proposed by Loudon [13], and applied by Caves 
and Crouch [14] to model distributed squeezing losses in a single-mode parametric amplifier. 
We present a fully analytical model of the relative squeezing produced by this system in the 
continuum limit, including losses on both beams. 
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Fig. 2. Competing gain and loss processes modelled by interleaved stages of squeezing 
(SQZ) and loss. 

Each stage comprises ideal squeezing (by parameter s) followed by loss (represented by 
transmission coefficients on the probe and on the conjugate). Consecutive stages of the 
model are related by combining Eq. (2) with Eq. (4) to give 



On+i = fa(coshi fl„ + siuhi bij) + \ / I — x„ 



1+1, 



K+1 = ffc(sinhi dn + coshi b^ + Jl-tl yl_^_^ 



(6) 



The overall model is parameterized in terms of the overall squeezing parameter S in the absence 
of losses, and transmissions Ta and Ti, in the absence of squeezing. These are related to the 
incremental coefficients above by s = S/N, = T^^'^'^ and f/, = T^^'^^ . 

This transformation can be written in matrix form and applied recursively to express the 
output beam operators qn and b^j as a sum of the incident beam operators (aq and Pq) and the 
injected vacuum operators (i,- and y\) as 




where 



fa cosh i fosinh,? 
ffosinhi ficoshi 



(7) 



(8) 



For algebraic simplicity, we introduce a polymorphic operator £, consisting of the probe anni- 
hilation and conjugate creation operators z\ = flo and Z2 = b^, followed by the injected vacuum 
operators zn+i = Xi and Z2i+2 = for 1 < ' < A^- Expanding Eq. (7) in terms of this operator 
and a set of coefficients a, and )3, gives 



(9) 



ON = aiflo + 0:2^0 + (^i'^l + 0!4i'| H (- 0:2^+1-% + CClN+iylj = Y.i 

K = Pi^o + Ihbl + jSsii + j34.yl H h Pin+ixn + PiN+iyl/ = A^'- 

Using standard statistical identities, the relative intensity operator and its variance are 

Na - Nh = aljON - b^bN = 1^,- - Pilij)z}zj + 1 and 

Var (Na-Nb) = ^ Aj3,)(a^a, - j3^j3,) CoVar (z]"!,-, z^z,) . 

In the bright probe beam approximation (^0 ^ I), this variance simplifies to 

Var (M, - Nh) = (a? - j32)2var {No) + ^ (ai a, - jSi A)' (^o> • 

The shot-noise limit is (A'a +Nb) ~ (oc^ + l5i){No), so the degree of squeezing for an initially 
shot-noise limited probe with Var (Nq) = (Nq) is 

2 



y^{Na-Nt) _ l.l,ia,at-m) .... 

It remains to express the a,, j3, coefficients in terms of the model parameters T^, Tf, and S, and 
hence obtain an ab-initio expression for the degree of squeezing. 
Equating the coefficients of Eq. (7) and Eq. (9) leads to 

f]=X« and h'^' 

Pi PlJ \P2i+l P2/+2 




Hence each of the vacuum modes i, contribute a term to the variance in Eq. (10): 
{a,a2i+i-mi+if - (ai, -pi) (^^':+') (a2,+i, Pl+i) "^J 

Each operator contributes a term analogous to Eq. (11), but whose diagonal matrix is 
[0 1 — f^J. Temporarily neglecting the ; = 1,2 contributions to the variance in Eq. (10), and 
summing over the vacuum contributions gives 




The continuum behaviour is recovered in the limit oo. To obtain a closed form expression 
for the sum, the infinitesimal parameters are expanded as a power series in 1 /N. Expanding the 
elements of A in Eq. (8) gives 

A.l + lA„ + 0(^) Where A„.(i">f" ^ J . ,13) 

Similarly taking tf^ = exp(^logra) ~ 1 + ^ logT,,, the sum in braces in Eq. (12) is 

X = V A^-' (It) A'^-' = - V A' T A' with T=('^J^" ° ^ V 

It can be easily verified that this sum obeys the geometric series relation 

AXA - X = -1^ (A^ T A'v - T). 

Expanding to order 1 /N using Eq. (13) gives 

AXA-X= if{AoX + XAo} => AqX + XAq = A'^T A'^-T. 

Taking the limit N ^ °°, the neglected 0(1 /N^) terms vanish and A'^ exp(Ao). Hence the 
sum X converges and obeys 

AoX + XAo = exp(Ao)T exp(Ao)-T. (14) 

This is a system of four linear equations for the elements of X in terms of the model parameters 
Ta, Th and S, and can be solved algebraically. 

The sum X contains all terms in the variance of Eq. (10) except / = 1,2 which correspond to 
the probe and conjugate coefficients. The probe contribution is 

while the conjugate contribution (ai a2 — PiPi)^ has diagonal [0 IJ . Computing the full vari- 
ance sum in Eq. (10) yields 



Predicted squeezing with probe absorption 
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Intrinsic gain, G = cosh^ S 

Fig. 3. Predicted squeezing for four-wave mixing with negligible conjugate absorption and 
net detection efficiency 77 = 85%. 



We introduce one final stage of loss to model optical losses after mixing (as in §4), scaling each 
coefficient by the relevant transmission factor {^/f\^ or ^/f\h) and introducing the extra vacuum 



contributions r\a{l — i7fl)oc^ + ~ lh)Pi- The net variance in Eq. (10) is therefore 



Var(^„-^i) = (ai, -/3i) {p(e2Ao+x)P+(l-P)p} (^0) where P= (^'J 



ria 
lb 



The measured beam powers relative to the incident probe power (the "effective gains") are 

and C..||.,.<« with Q^^*"®- (15) 

and the relevant shot-noise Hmit is Var {Na-Nh)^^^^ = {N,,) + {N/,) = (r;„af + rihj5f){NQ). 

Evaluating all the contributions to Eq. (10), we obtain an analytic expression for the degree 
of relative intensity squeezing produced by this system. This expression is algebraic in r„, Tj, 
and S, but runs to a dozen typeset lines. However, special cases are readily derived and provide 
physical insight not readily accessible from numerical models. 

In the experimentally studied case (Refs. [1-7]), detection efficiencies are carefully balanced 
("Ha = 'nb = 'n) ^nd the far-detuned conjugate experiences negligible absorption (T/, = 1). The 
corresponding degree of squeezing is 

25sinh2£ _ S \og^T„smh'^t 

'<^cosh 2^+;t "2^3cosh ' 



with parameters ^ = i ^ 1 65^ + (log r„ - log 7), ) 2 and tanhj = (log r„ - log 7},)/4^ . The three 
terms describe the shot-noise limit, coiTelations from four-wave mixing, and injected vacuum 
noise. 

Figure 3 shows the noise figure as a function of the probe transmission Ta and intrinsic 
mixing gain G ~ cosh^^ (in excellent agreement with the numerical model of Ref. [3]). Note 
the counter-intuitive result that the strongest squeezing is obtained with imperfect transmission 
(Ta < 1). This is because the shot noise carried by the incident probe beam is also amplified 
by the mixing process, and a small "optimal" level of probe loss decreases this contribution to 
the measured noise power before injected vacuum noise dominates. This optimal level is easily 
obtained by minimizing Eq. (16). 



A. "Reverse" double-A mixing configuration 



B. Predicted squeezing for losses on conjugate 
and probe 
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Fig. 4. (A) Four-wave mixing in tlie "reverse" configuration. (B) Predicted squeezing for 
four- wave mixing in the normal/reverse configurations with intrinsic gain G = 3. 

The model can also be applied to other four-wave mixing systems. For example, interchang- 
ing the probe and conjugate wavelengths produces the "reverse" configuration (Fig. 4A). This 
system is interesting as the weakly-coupled conjugate transition is brought closer to resonance, 
producing much higher intrinsic mixing gain for the same beam powers. In this case, the probe 
experiences negligible absorption {Ta = 1) compared to the conjugate (T^, < 1), and the pre- 
dicted squeezing is 



The vacuum noise term in Eq. (17) is considerably larger than in Eq. (16), resulting in several 
decibels difference for moderate levels of absorption (Fig. 4B). Unlike the probe losses dis- 
cussed above, losses on the conjugate only destroy correlations and introduce noise, so squeez- 
ing by four-wave mixing in the reverse configuration is always less effective for the same level 
of intrinsic mixing gain. 

Finally, it is worth noting for consistency that in the limit of both r„, ^ 1, the post-mixing 
optical-loss result of Eq. (5) is obtained. 

5. Experimental diagnostic 

The analytic model derived above provides a simple yet powerful tool for optimizing the de- 
gree of squeezing obtained experimentally. It not only predicts which parameters will provide 
optimal results, but can also be used as a diagnostic tool to determine which factors are limiting 
the experimentally measured degree of squeezing. 

To demonstrate this, we constructed a four-wave mixing apparatus as described in Ref. [1]. 
A 400 mW pump beam intersects a 100 /iW probe beam with 1 /e^ beam waists of 630 /im and 
375 jim respectively at an angle of 0.3° within a pure ^^Rb vapor cell of internal length 7 mm 
heated to 130°C (Fig. IB). The probe was generated by an AOM with fixed detuning 3040 MHz 
below the pump, which was scanned across the Doppler broadened D\ resonance at 795 nm. 
The relative intensity between probe and conjugate was measured with a balanced photode- 
tector (Thorlabs PDB150A), refitted with high-efficiency photodiodes (Hamamatsu S3883, net 
efficiency 95%). The overall detection efficiency of the system was 77 = 85 ± 1%. The relative 
intensity noise was measured with a Rhode & Schwarz FSP7 spectrum analyzer at an analysis 
frequency of 1 MHz with 30 kHz resolution bandwidth. 

The effective probe and conjugate gains (Gq and G^) were measured as a function of pump 




(17) 



A. Measured Probe and Conjugate Gain C. Measured and Predicted Squeezing 
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B. Computed Gain and Transmission Parameters D. Measured and Predicted Noise Power 
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Fig. 5. (A) Measured probe and conjugate gain across four-wave mixing resonance and 
(B) inferred model parameters. (C) Experimentally measured squeezing and (D) associated 
noise powers compared to model predictions for these parameters. Detunings are for pump 
beam and are measured above the centre of the ^^Rb 55i/2(^ = 2) 5Piy2 transition. 

beam detuning (Fig. 5A) and used to simultaneously solve Eq. (15) for the intrinsic gain G and 
probe transmission Ta (Fig. 5B) via the coefficients ai and a2- Note that the gain resonance 
extends well into the Doppler-broadened absorption resonance for detunings below 600 MHz 
from the line-centre, demonstrating strong competition between the processes. 

The model predicts significantly stronger squeezing should be possible for these parameters 
than measured experimentally (Fig. 5C). As our incident probe beam was measured to be shot- 
noise limited, this demonstrates that our noise measurement was limited by technical factors 
and not by insufficient intrinsic gain or excessive absorption. 

One cause for our discrepancy is the technical difficulty in eliminating the bright pump beam 
after the vapor cell. Detection of the pump introduces uncorrelated fluctuations and increases 
the measured noise level [3]. This is evident around zero detuning, where the probe is fully 
absorbed (Fig. 5A) but the measured noise is far above the standard quantum limit (Fig. 5D). 
Upon blocking the pump beam, the relative intensity noise was measured at the shot-noise 
limit. Subtracting this background level implies that —4. 1 dB of squeezing would be obtained 
by eliminating cross-beam detection. 

The model assumes that gain occurs uniformly throughout the vapor cell, which requires that 
the beams be overlapped over the entire region. If the beams are not properly overlapped, the 
mixing strength decreases and gain becomes spatially varying. As losses are unchanged, the 
relative intensity noise increases as a result. This is likely the cause of the remaining discrep- 
ancy, and improved squeezing could be achieved by manipulating beam alignment and waists 
to ensure proper overlap throughout the cell. 



6. Conclusions 



We have presented a method for analytically calculating the degree of squeezing produced by a 
four-wave mixing system in the presence of absorption. Our model included the contributions 
from arbitrarily many injected vacuum modes that were subsequently squeezed by the system, 
producing an ab-initio quantum mechanical description of the introduced losses. Our general 
result is not reliant on implementation details and can be applied to analyze any four-wave 
mixing scheme, while the matrix methods techniques we developed can be applied to model 
other systems both within quantum optics and more generally. 

We considered two special cases, corresponding to the experimentally studied system (in- 
troduced in Ref. [1]) and the "reverse" configuration, with probe and conjugate wavelengths 
interchanged. We presented closed-form expressions for the relative intensity noise in these 
cases, and demonstrated that a small level of probe loss was desirable to suppress amplification 
of the initial shot-noise. The reverse configuration was shown to produce the same squeezing 
at equal gain for ideal transmission, but was significantly more sensitive to losses. 

The model was applied to analyze experimentally measured squeezing and determine the 
intrinsic mixing gain and transmission factors of the four-wave mixing resonance. Compar- 
ing the expected squeezing to measured results provided insight into the factors limiting our 
measurement and hence where to direct effort in optimizing the many free parameters of the 
system. 

While it should be noted that this model considers the propagation of a Gaussian beam mode 
only, an arbitrary beam can be analyzed as a product state of orthogonal spatial modes, with 
each mode independently squeezed [4]. The model established in this paper may therefore be 
applied to each pair of spatial modes and the resulting noise powers summed to obtain the 
overall relative intensity squeezing for a multi-mode beam. Such multi-mode squeezing has 
been experimentally demonstrated [5], with promising applications in quantum imaging. 



